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SUMMARY 


This paper presents the algorithms and logic for use in the implementation 
of instrument landing system (ILS) localizer signals for the generation of 
navigation information for initial approach guidance. Navigation position 
estimates are based on position updates determined from range informa- 
tion measured from a DME station (not necessarily colocated v/ith the ILS 
antennas) and ILS localizer angular information. ILS volumetric coverage 
and DME geometric checks are performed to ensure that proper radio navigation 
inputs are being utilized. These algorithms and volumetric and geometric 
checks are designed so that they could be added to area navigation systems 
utilizing position difference information to determine navigation estimates 
with minimum software modification. 

INTRODUCTION 

Flight tests with the NASA TCV B-737 airplane (ref. have shown that 
its area navigation system can provide new guidance information presented 
on cathode ray tube (CRT) displays that results in improved lateral and 
vertical flight path tracking during an approach to landing. Additional 
operational experience (ref. 2) has shown that the area navigation system 
can also provide guidance for flight along curved paths to avoid noise 
sensitive areas during an approach to landing. To complement the approach 
to landing guidance, the navigation computer must provide an accurate 
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estimate of the airplane's position. 

The TCV area navigation system provides an estimate of position from 
position updates derived from various radio navigation facilities and 
inertial navigation system derived velocity. It normally uses dual DME radio 
inputs for enroute navigation and has used microwave landing system signals 
for precision navigation to provide approach and landing guidance. To ensure 
that the best possible position estimate is obtained while operating in an 
area of instrument landing system ( I L S ) signal coverage, algorithms and logic 
have been developed that will utilize localizer signals and range information 
from a single DME station. 

The purpose of this report is to describe the equations and logic used 
to generate a navigation position estimate in the ILS localizer signal 
environment. The equations used to calculate position difference components 
derived from 1) ILS localizer deviation measurements and randomly located 
DME range measurements and, 2) only ILS localizer deviation measurements are 
described. A summary of the NASA TCV B-737 airplane's estimate process will 
also be described so that the total navigation process may be shown. 

SYMBOLS 

A - distance between the DME and ILS localizer antenna - n.mi. 

a - East component of A - n.mi. 

b - North component of A - n.mi. 

D - DME corrected for slant range - n.mi. 

O' - DME slant range - n.mi. 

DME - distance measuring equipment 
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- position difference vector 

- component of the position difference vector, DP, perpendicular 

to the runway centerline 

- magnitude of vector DPp - n.mi. 

- cycle interaction time of navigation computer - sec 

- ellipticity constant, 0.003367 

- altitude of aircraft above mean sea level - m 

- altitude of DME antenna above mean sea level - m 

- altitude of ILS localizer antenna above mean sea level - m 

- unit coordinate vectors 

- navigation update mode; inertial velocity, ILS localizer, and 

DME 

- instrument landing system 

- navigation update mode; inertial velocity and ILS localizer 

- position update gain 

- velocity update gain 

- outer radial limit of localizer volumetric check - n.mi. 

- inner radial limit of localizer volumetric check - m 

- axes of orthogonal coordinate system oriented towards True 

North 

- angle formed by the vector Z and a line between the airplane 

and DME antenna - deg 

- Earth radius - m 

- Meridional radius of curvature - m 

- Normal radius of curvature - m 

- unit vector perpendicular to the runway centerline 
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North and East components of system velocity estimate - kts 

North and East components of inertial ground speed - kts 

axes of orthogonal coordinate system oriented along the 
runway centerline 

coordinates of vector Z transformed into the X', Y' 
coordinate system 

airplane's estimated position vector 

magnitude of Z the airplane's estimated position vector - 
n.mi. e 

North and East component of 1 the airplane's estimated 
position vector - n.mi. 

airplane's measured position vector 

magnitude of Z” M - n.mi. 

North and East components of Z M 

vector of airplane's position estimated radially along the 
measured localizer deviation angle 

magnitude of position vector Z^ - n.mi. 

North and East components of position vector Z^ 

angle formed by the DME antenna, origin, and measured 
airplane position - deg 

angular limit about angle P where DME signal usage will be 
inhibited - deg 

North and East components of position difference -n.mi. 

North and East components of system velocity update - kts 
latitude and longitude update estimate - deg 
localizer deviation angle - deg 
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p - relative angle between the DME & ILS localizer antenna - deg 

^DME ’ *DME " ^ at i tuc ^ e » longitude of DME antenna location - deg 

<f> e , - latitude and longitude of airplane position estimate - deg 

<{> 0> - latitude, longitude of ILS localizer antenna location - deg 

ij; - runway heading relative to true north - deg 

q - vertical angular limit of ILS localizer volumetric check - deg 

a) - lateral angular limit of ILS localizer volumetric check - deg 

Subscripts: 
t - time 

DISCUSSION 

NAVIGATION POSITION ESTIMATE DESCRIPTION 
Position Difference 

The TCV B-737 navigation computer is software controlled (ref. 3 and 4) 
to select and tune two appropriate DME and/or VOR stations in the vicinity 
of the airplane. A position difference between the point defined by the distance 
and/or azimuth information received from the stations and the previous position 
estimate is determined. This position difference is divided into North and East 
components to be used to calculate a new position estimate. 

Other sources of navigation information may also be used to determine 
position differences. While operating in the ILS localizer environment, these 

position difference components are calculated with localizer and DME information. 
This estimation mode is shown on the pilot's electronic map display as ILD 
(inertial, localizer, DME). If no DME is available, the estimation mode is ILX. 
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Position Estimate 


In the TCV B-737 airplane's navigation system, the same position 
estimation calculations are used regardless of the manner in which the 
position difference is determined. The first step in the position estimate 
process is to develop a system velocity update from the position difference. 

""t - ♦ K 2 

AV = AV r + K 9 AP f 

4 h t-i 4 


A system velocity estimate is obtained by summing the system velocity 
update with ground speed obtained from the inertial navigation system. 





V c = AVp + V F 

4 t 4 


A position update, in terms of latitude and longitude, is obtained using 
the system velocity and the position difference as follows: 

A* t = (V„ t At + K, AP Nt >/R„ t 

AX t = (V £ At + K, AP e )/R n . 

t t t 

This position update is based on an oblate spheroid earth model (ref. 3) by 
using appropriate radii of curvature in the North/South and East/West 
directions. 

R M = H a/c + V 1 - 2F + 3F sin2 4> e ) (North/South) 

r m = [H. /c + R E (1 + F sin 2 4 )] cos <fc (East/West) 

t a/ 4 b e t-i 4-1 

The updated position estimate is found by summing the previous 
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position estimate with the position update terms. 


<P e ~ <P e 

e t e t-l 
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ILS/DME (ILD) POSITION ESTIMATE 
General Solution 

Figure 1 shows the geometry of the ILS antenna, the DME antenna, the 

airplane's position estimate and its position as measured with ILS localizer 

and DME information, and the runway with an extended centerline. An 

orthogonal coordinate system, with its origin placed on the ILS localizer 

antenna, is oriented in a True North direction. Since the DME may be either 

automatically or manually selected within geometric constraints relative to the 

airplane, it is not required to be co-located with the ILS localizer antenna. 

By determining the relative geometry between the ILS localizer antenna 

and the DME and the airplane's position measured with the localizer and DME 

information, the vector between the origin and the airplane's measured 

position, Z^, may be found. The vector between the origin and the airplane's 

previously estimated position, Z » may be found directly from their known 

e t-l 

latitudes and longitudes. The desired position difference used to form a new 

position estimate may then be found by subtracting from Z g 

t t“l 

Calculation of the Airplane Position Vector T M 

t 

Figure 1 shows that Z M is one side of the triangle formed between the 

l t 

origin (ILS localizer antenna), the airplane's measured position, and the 
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DME antenna. Known quantities used to determine 7^ in this triangle include 

the runway heading, ^ R ; the slant range of the DME, D'^ the ILS localizer 

deviation angle relative to the runway centerline, n t ; and the latitudes 

and longitudes of the DME antenna, <f> DME , A DME ; and the ILS localizer antenna 

location d> . A . 

o o 

Figure 2 shows the angular geometry and distance between the origin 
and DME antenna, A. The length A and the relative angle between the DME 
antenna from the North axis, y, remains constant and is calculated only 
once. If a new DME is tuned, then A and y are recalculated. 

The length of A is determined by vectorially summing its components 
a and b. 


a = (Xce - A o )(60) c.s *L^SSL 


b = ^ DME ' * O X 60 > 


-1 


a 2 + b 2 


The angle y is found by 


y = tan 


-1 


7 » 


o < y < 2 7T 


n.mi . 

n.mi . 
n.mi . 


The angle ou, formed by side A and vector Z„ 
movement and should be continuously recalculated, 
between o and tt, inclusive, and is found by 


a t = 


2 tt -V + (v'p - 


, may change with airplane 
Angle a^. has a range 
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If the absolute value of c^is greater than tt , then: 


cl 


2tt - ol 


The magnitude of the DME reading as measured in the airplane, D'^, is 
the slant range distance between the ground-based DME antenna and the airplane. 
This distance was slant-range-corrected to determine the ground distance, D t , 
between the airplane and DME 


Dj. = D't sin 


cos 


-1 


H a/c t " H DME 


D' 




Angle Fj. is formed by the vector Z^and side D^. This angle may 

continuously vary and is calculated 20 times per second. Knowing the relation 

D 4 
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angle P t i s determined by 

r A 
~ D 


P t = sin -1 


i not J 


It must be determined if angle P t i s an obtuse or an acute angle. This 
is accomplished by comparing the square of side A with the sum of the square 
of side Dj. and the square of the magnitude of the estimated airplane position 
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The magnitude of Z is used as an approximation for Z M . The magnitude 
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. “ 0 ( 60 ) 
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ZeE t-i = ( S-, ' x ° )(60) cos *° * Vi 


The magnitude of Z M is found by 

”t 

Z^= Aj-cos (c^)+ D^cos (P t ). 

The North and East components of Z^are found knowing the angle, 
toR-nj-), between the North axis and Z M . 

Z M t C0S %-”t> 

^ Z (^. sin ^R -r V 

(Calculated of the Position Difference in North and East Components: AP 

AP e -ILD Update Mode 

The position difference, in North and East components, is now found 
by subtracting the North component of Z from the East component of Z„ and 
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t-1 


"t 


the East 


component of Z g ^ from the East component of Z^ . 
- Z. 


A V H-i 


-MN 4 


AP E t " ZeE t-l " Z ME t 

These position difference components are then used directly in the 
navigation position estimation algorithms described previously. 


ILS LOCALIZER ONLY (ILX) POSITION ESTIMATE 
General Solution^ 

In the event that a DME can not be tuned, the navigation computer will 

utilize the ILS localizer signal and inertial velocity to determine a position 

Vhe ILX position estimation mode was developed by the Boeing Commercial Airplane 
Company at the NASA's request. 



estimate. Figure 3 shows the geometry of the ILS localizer antenna, the 
airplane's position estimate and its estimated position on a measured localizer 
deviation angle, and the runway with an extended centerline. An orthogonal 
coordinate system, with its origin placed on the localizer antenna, is 
oriented in a true north direction. 

Radio position difference in the ILX update mode is limited to a direction 

perpendicular to the runway centerline. Since no DME information is 

available, a radio position difference cannot be developed in a radial direction 

from the localizer antenna. However, inertial velocity is still utilized 

in the radial direction and will supply inputs for a new position estimate 

in the position estimate algorithms. 

North and East position estimates are found in the following manner. 

A position difference vector, Dfy, is found by subtracting the estimated 

position vector, Z , from the estimated position vector on the measured 
e t-l _ _ 

localizer deviation angle, Z D • The component of DP perpendicular to the 

K t 1 

runway is found and broken into North and East components. These components 
are used in the navigation position estimate algorithms. 


Calculation of DP^. 

To determine it is necessary to calculate the position estimate 


vector, Z g ^ ^ , in North and East components from the latitudes and lonaitudes. 
of the last position estimate and the ILS localizer antenna location. 
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A vector Z^,of the airplane's estimated position along a measured 

localizer deviation angle is determined in the following manner. Since no 

radio updates can be obtained in a radial direction from the localizer 

antenna, it will be assumed that the estimated radial distance from the origin 

is correct. Hence, the vector length of Z D and Z will be the same 

_ R t e t-l 

(directions may differ to obtain DP t ). 

is: 


The length of Z D and 2 

K t e t-i 


\ s Vf 




■1 


t-1 


The North and East components of are found knowing the angle, 
(^R -r t ) between the vector and the North axis. 


= Z 


R t C0S (W 
‘RE. = Z R. Sin ^R"V 

t C A 
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DPMI's found vectorially 


DP*. - l n - 1. 


c t-l 


Calculation of the Position Difference: AP^ , AP^ - ILX Update Mode 

t t 

The magnitude of the component of perpendicular to the runway 
centerline is obtained by the vector dot product of D^with a unit vector, 
u, perpendicular to the runway centerline. This results in the magnitude 
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of the position difference. This magnitude is then multiplied times 
the unit vector to obtain the position difference in North and East 
components. These components are then used directly in the navigation 
position estimate algorithms. 

A 

The unit vector, u, shown in figure 3, is: 

A A A 

u = - cos (i^ R ) i + sin (ij/ R ) j 

The magnitude of the component of D(£ perpendicular to the runway 
centerline is: 

As 

0p p t = op t ' u 

DP p t ‘ ' " Z eE t _ 1 ) cos + ^ Z RN t " Z eN t _ 1 * s,n **R* 

The North and East components of the position difference are: 

A \ = DP p t sin <*R> 

AP^ - -DPp cos . 

These position difference components are then used directly in the 
navigation position estimate algorithms. 

ILS LOCALIZER/DME/AIRPLANE POSITION VALIDITY CHECK 

Since random DME selection, automatic frequency tunning, and other 
means of automatic software control will be utilized in the navigation 
computer, a check must be made to ensure that the navigation computer is 
using appropriate navigation data that will result in accurate position 
updating. Obviously, if improper navigation data are being utilized, position 
estimates based on that data must be inhibited. 
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To preclude the possibility of using the wrong ILS localizer navigation 
data, a volumetric geometry check is made to determine if the airplane's 
position estimate is within the localizer azimuth boundaries. If the 
airplane is not within these boundaries then ILS updating is inhibited. 

These boundaries include lateral and radial limits of coverage as shown in 
Figure 4. The lateral limit of the localizer angle of coverage is +to from 
the back azimuth antenna. Radial limits require that the airplane be within 
L n.mi. of the localizer antenna, but not closer than M meters. A vertical 
angular limit of n degrees from the localizer antenna is also applied to 
ensure that the airplane's altitude does not exceed the localizer 
signals vertical coverage. Localizer back course information is not used 
for updating since all back course signals will not provide valid navigation 
information. 

A new orthogonal coordinate system, with its origin located at the ILS 

localizer antenna and X' axis parallel to the runway centerline, was used to 

make the geometric checks. To make the checks, the North and East components 

of the estimated position vector 1 , are transformed into the new X', Y' 

e t-l 

coordinate system as follows: 


X 

Y' 


' = Z « cos (i|> R ) + Z E 

e t eN t-l K et t-l 


p Z eN sin ' Z eF 

e t eN t-l K et t-l 

The lateral azimuth check is then 


-X' tan .. < Y' < X' tan to. 
^ “ e t“ h 

The radial check is 

M meters < X' < L n.mi . 

- e t“ 

The vertical check is 

“a/c^ H o + X V an 


sin (t R ) 
cos (<jt R ). 
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A geometric check (figure 5) of the DME ground station's location relative 
to the airplane and runway centerline is continuously performed to prevent 
the possibility of large position estimate errors caused by the DME inputs. 

As the angle P t formed by the localizer antenna, airplane, and DME ground 
station, approaches 90°, small errors from the DME range information can 
result in significantly larger radial position errors along the localizer 
deviation angle on which the airplane is located. Hence, if angle falls 
within a prescribed angular boundary, 3, about P=90°, the DME signal will 
not be used and a new DME will be tuned or ILX updating will occur. This 
geometric check is defined as 

(§ + 6) .< P t < (| - B). 


CONCLUDING REMARKS 

Flight operations of the NASA TCV B-737 airplane have shown the benefits 
of using the area navigation computer to provide guidance to fly within the 
ILS area of coverage. Its area navigation system may use various radio 
■inputs to calculate a position estimate for course guidance. 

ILS localizer deviation signals may be utilized on the TCV B-737 airplane 
by algorithms that generate position difference information for i/se in the area 
navigation system. These algorithms are designed so that they may be 
added to area navigation systems utilizing position difference information to 
determine navigation estimates with minimum software modification. 
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Izer, DME, airplane estimated and measured positions, and runway with 
led centerline geometry— ILD update mode. 





No DME 
update 



Figure 5.- DME ground station geometry check for ILD position updating 



2. Government Accession No. 


3. Recipient's Catalog No. 


1. Report No. 

NASA TM-80167 


15. Supplementary Notes 


16. Abstract 

This paper presents the algorithms and logic for use in the implementation of 
instrument landing system (ILS) localizer deviation signals for the generation 
of navigation and guidance information. Navigation position estimates are 
based on range information from a randomly chosen DME and ILS localizer deviation 
information. ILS volumetric coverage and DME geometric checks are performed 
to ensure that proper radio navigation inputs are being utilized. These algorithms 
and volumetric and geometric checks were designed so that they could be added to 
most existing area navigation systems with minimum software modification. 


’ For sale by the National Technical Information Service. Springfield. Virginia 22161 


18. Distribution Statement 

Unclassified - Unlimited 

Subject Category 04 


21. No. of Pages 

22. Price* 

21 

$4.00 


17. Key Words (Suggested by Author(s)) 

navigation 

instrument landing system 
area navigation 


19. Security Classif. (of this report) 20. Security Classif. (of this 

Unclassified Unclassified 


4. Title and Subtitle 

ALGORITHMS AND LOGIC FOR INCORPORATING ILS LOCALIZER 
INFORMATION INTO THE NASA TCV B-737 AIRPLANE AREA 

NAVIGATION SYSTEM 


7. Author(s) 

Charles E. Knox 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


5. Report Date 

September 1979 

6. Performing Organization Code 

8. Performing Organization Report No. 

10. Work Unit No. 

513-52-13-52 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 

14. Sponsoring Agency Code 

















